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pH ■ Abstract 

'-^ , Recently, a new type of battery has been proposed that rehes on the principle 

^ I of self-assembling of a liquid metalloid positive electrode, a liquid electrolyte, 

^ I and a liquid metal negative electrode. While this configuration has been 

c/j ■ claimed to allow arbitrary up-scaling, there is a size limitation of such a sys- 
tem due to a current- driven kink-type instability that is known as the Tayler 

^! instability. We characterize this instability in large-scale self-assembled liquid 

^ ' metal batteries and discuss various technical means how it can be avoided. 
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^ . 1. Introduction 



o 



With growing deployment of intermittent renewable energy sources, such 
as wind and solar, large scale electricity storage becomes an issue of increasing 



in 

O ' importance. 

^ ■ In a recent patent application [ij (see also |2| for more details) a new 

type of battery has been proposed that has the potential to become a key 
ingredient in balancing supply and demand of electrical energy. Its working 

r> I principle relies on reversible ambipolar electrolysis within a self-assembled 

c^ ; layered structure of a liquid metal, an electrolyte, and a metalloid. Histori- 

cally it is worth to note that the idea of self-assemblage due to appropriate 
density differences was already discussed in |3|. Such an assembly of three 
purely liquid layers, without the interference of any solid phase, allows for 
the highest possible reaction rates known in electrochemistry. The absence 
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of any vulnerable ceramic electrolyte, as it is used in sodium-sulfur (NaS) 
batteries j^, makes the liquid metal battery a promising candidate for al- 
most unlimited upward scalability, which appears to be essential for economic 
competitiveness. 

In this paper we will discuss a possible limitation of the upward scalabil- 
ity of liquid metal batteries due to magnetohydro dynamic (MHD) instabil- 
ities in conducting fluids under the influence of externally applied electrical 
currents. The first problem that comes to mind here is the interfacial insta- 
bility that is well known in aluminium reduction cells p, l6| . This interfacial 
instability starts from small interface deformations between cryolite and alu- 
minium which lead to a mainly horizontal disturbance current in the liquid 
aluminium. The Lorentz force resulting from this horizontal current in in- 
teraction with a given vertical magnetic background field (an unwanted side 
effect of the bus bars supplying the current to the cell) drives a horizontal 
motion with a rotating interface that under certain conditions may become 
unstable. 

The onset of this instability depends on the height of the cryolite layer 
which must not fall under a certain value. The usual way to avoid this 
instability is, therefore, to use a thicker electrolyte layer which results in a 
higher electrical resistance. Actually, this resistance is responsible for the 
giant Joule losses in aluminium reduction cells which are responsible for the 
consumption of about 2 per cent of the electricity generated world-wide |6|. 
Indeed, this interfacial instability could play a role in liquid batteries, too, 
and it should be considered carefully in order to determine a lower threshold 
for the thickness of the electrolyte layer. 

The point of the present paper, however, is to highlight and character- 
ize another type of instability which could set a serious limit for the upper 
size of liquid metal batteries. This instability is well known in astrophysics 
under the label Tayler instability (TI) |7| (sometimes also called Vandakurov- 
Tayler instability |8|). The TI can be considered as a limiting case of the 
well known kink-instability in plasma physics that occurs if the so-called 
Kruskal-Shafranov stability criterion is violated, i.e. when the ratio of axial 
to azimuthal magnetic field falls below some critical value |9|. 

In our context the TI is a kink-type (i.e. non-axisymmetric) instability 
that occurs if the current through a column of a liquid metal exceeds some 
critical value in the order of kA, depending on the combination of the ma- 



terial parameters conductivity, viscosity and density [10|, lUl ■ If this current 



threshold is exceeded, the TI would lead to a quite vigorous motion and very 



likely to an undesired mixing of the three layers which should definitely be 
avoided for the battery to function properly 

The main goal of this paper is to figure out the principle importance of the 
TI for liquid metal batteries. In a very first attempt to quantify the critical 
current, we will use a one-dimensional equation system for the hypothetical 
case of a cylindrical fiuid with infinite length. This is of course a strong 
idealization of a real battery with its finite ratio of length to radius and 
its possible deviation from the cylindrical shape. Note, however, that there 
are indeed visions of high capacity batteries for which the ratio of length to 
radius of the liquid metal and/or metalloids might be in the order of one or 
larger [l|, |2|. For any detailed battery design the correct determination of 
the current threshold would require expensive two-dimensional simulations 
(in the spirit of [12]) which is beyond the scope of the present paper. 

The paper starts with a short outline of the recent ideas about liquid metal 
batteries. It continues with a characterization of the TI in ideal fiuids which 
will help us to find a most simple means on how the TI can be avoided. Then, 
for the most widely discussed hquid metals Mg and Na, we will compute the 
critical currents which should be considered as conservative estimates for the 
onset of the TI. For those liquids we will also propose a simple means by 
which the TI can be shifted to significantly higher critical currents. 

2. Liquid metal batteries 

Liquid metal batteries are thought to play a big role in future grid energy 
storage. Sodium-sulfur batteries with a power around 30 MW are already 
used in a Japanese wind-park. They work with a cell voltage of around 2 V, 
at temperatures of around 300°C, and they have AC- AC efficiency of 75-80% 
[^. Their main problem, however, is the use of a ceramic /3"- alumina solid 
electrolyte which is susceptible to thermal shocks and which also results in 
size limitations due to the difficulties to manufacture it in thin cross sections 
while maintaining structural strength. 

In order to avoid problems of this sort, a perfectly liquid battery was 
proposed in [l|, |2[. Such a battery consists basically of three liquid layers 
which are self- assembling due to their different densities and their mutual 
immiscibility (see Fig. [T]). The lowest layer, consisting of a metalloid like 
Sb or Bi, represents the positive electrode of the battery where the negative 
ions of the metalloid are oxidized during the charge process. The second 
layer, the electrolyte, is a molten salt which is chosen in such a way that it 



is immiscible with both the metal and the metalloid, allows for the diffusion 
of their respective ions, is electronically insulating, and ionize the molten 
product of the metal and the metalloid. As potential electrolyte components 
Na2S, Li2S, CaS, Na2Se, Li2Se, and NaaSb are discussed J2|. Obviously, the 
electrolytes' conductivity should be as high as possible in order to minimize 



Joule losses. However, unlike in liquid metal concentration cells [13[, where 
the electrolytes' role is merely that of an ion conductor, liquid metal battery 
electrolytes have to absorb a considerable amount of metal and metalloid ions. 
This requires high solubility of the metal and metalloid ions in the electrolyte 
as well as a relatively large electrolyte volume. The latter exigence entails 
thick electrolyte layers opposing the need for low cell resistance. The use of 
external electrolyte reservoirs and a circulation system would be a possible 
mean to meet these conflicting goals. 

The third and top layer consists of a molten earth alkaline metal such as 
Mg or, possibly, an alkaline metal such as Na. This layer serves as negative 
electrode where the metal cations are reduced during the charging process. 

It goes without saying that the safe operation of such a cell requires the 
self-assembled layer structure of the three liquids to remain stable during 
charging and discharging. 

For the quite similar case of aluminium reduction cells the interfacial 
instability is well known to set some lower limit for the thickness of the 
electrolyte below which the fluid starts to undergo a wavy motion that can 
ultimately interrupt the electrolysis. The same sort of instability has to be 
carefully considered for liquid batteries, too. 

The focus of this paper is, however, on the existence of an upper limit of 
the size of the battery beyond which the Tayler instability would stir up the 
liquids. 

3. The Tayler instability for the ideal fluid 

The Tayler instability is well known in astrophysics where it had been 
first discussed by Vandakurov |8| and Tayler [TJ. It can be considered as a 
limiting case of the kink-instability in plasma physics that sets in when the 
safety factor (basically the ratio of axial to azimuthal magnetic field) falls 
below some critical value. 

For an ideal fluid, without any viscosity and resistivity, Tayler |7| had 
given the condition of stability of a purely azimuthal magnetic field B^{r) 
with respect to non-axialsymmetric perturbations (proportional to exp (i(p)) 
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Figure 1: Conceptual design of a liquid metal self- assembling cell with coordinate system. 



in the following form: 



|-(r55(r)) > . 



(1) 



Despite its restriction to ideal fluids, which will be relaxed in the next 
section, it is worthwhile to analyze this expression a bit further. With view 
on our goal to avoid the TI in batteries, we will consider not only the obvi- 
ous case of a completely filled cylindrical vessel with a homogeneous current 
distribution, but also the more general situation of a homogeneous current 
through the liquid in a hollow cylinder, together with an additional indepen- 
dent current along the axis. 

Assume the fluid to fill the hollow cylinder between the inner radius r^ 
and the outer radius r^, with the ratio of the radii denoted hj rj = r[/ro. 
The current density j through this cylindrical gap is assumed homogeneous 
so that the total current through the liquid is If = n^rl — rf)j. Assume 



further the existence of an independent current I^ along the central axis of 
the cylinder. Applied to liquid metal batteries, this scheme is depicted in 
Fig. [21 Under these conditions the azimuthal magnetic field in the fluid 




Figure 2: Sketch of a liquid metal battery with a central bore and a conductor for a 
stabilizing axial current. 



acquires the form 



B^{r) = ar + (5/r 



(2) 



with the constants a. and /3 given by 

2ti r1 1 — rj 






^ 2. 






X — r^2 
The stabihty condition ([1]) can then be re-written in the form 



(4) 



/?' 



2 



.2 



^ - 2a/3 - 3aV^ > . (5) 

From the specific radial dependence of this expression it is evident that 
the most dangerous position for the instabihty to occur is close to Tq. There- 
fore, in order to identify the critical values for /3, it is sufficient to solve the 
quadratic equation 

13^ - 2aPrl - 3a V^ = . (6) 

only at Tq. This quadratic equation has two solutions which are Pi = —ar"^ 
and (32 = Sar^. Expressed in terms of the current /a and If, this can be 
rewritten into the following two conditions for stability: 

4 < -/f (7) 

or 



4 > /r^f^ . (8) 



3 + V^ 

These two inequalities are the main result of this section. They suggest 
two simple possibilities to avoid the TI. Either one applies 4 in opposite 
direction to /f: then 4 has to be at least as strong as If. Or one applies 4 in 
the same direction as If. Then 4 has to be factor (3 + ?7^)/(l — f]'^) stronger 
than If. 

This is illustrated in Fig. |3]where we show, for the particular case t] = 0.5, 
the (normalized) function —rB'^{r) for various values of /a//f- For all curves 
with —1 < la/If < 13/3 we can identify regions (at least close to r^) where 
the function —rB^{r) has an outward decreasing part. Both for /a//f < —1 
and for /a//f > 13/3, however, the function —rB'^{r) increases monotonically. 

For the purpose of batteries, the first version with /a < —If seems much 
more convenient. The limiting case /a = —If can be implemented by simply 
returning the battery current through the inner hole of the battery. 
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Figure 3: Dependence of the normalized function —rB'^{r) on the parameter /a//f for the 
particular case rj — 0.5. Note the monotonic increase of this function both for /a//f < —1 
and for la/If > 13/3 which indicates the absence of the Tayler instability. The thick 
dashed line with an arrow indicates the change of the profiles when going from negative 
to positive la/If 



4. The Tayler instability for real fluids 

While in case of an ideal fluid with a homogeneous current distribution 
the TI would set in already for arbitrary small currents, the main eflect of 
finite viscosity and resistivity is to shift this threshold to some finite value. 

In this section we will determine the critical currents for some real liquid 
metals which seem most important in our context, and we will study their 
dependence on the radii ratio t]. 

In the appendix we delineate the numerical scheme for the determination 



of the critical current which relies on the scheme outlined in |10l . Ill| . In 
this respect it is important to note that the relevant number here is the so- 
called Hartmann number Ha = BR\Ja/{pv)^ where B is the magnetic 
field strength and R is the radius of the cell. Hence, the larger the factor 
■\/o"/(pz/), the easier the TI will be excited for a given current. This means. 
in turn, that if we have a self-assembled layer of liquids, the most critical one 
of them is that with the highest factor. It turns out quickly that the typical 
metals (Mg, Na) forming the negative electrode are more prone to TI than 
the half-metals (Sb, Bi) forming the positive electrode. 

In the following we will consider two of the most relevant materials. The 
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first is Mg at 700°C with the parameters a = 3.61 x 10^ S/m, p = 1.54 x 
lO^kg/m^, u = 6.74 x lO^^mVs, the second is Na at 300°C with a = 5.99 x 
lO^S/m, p = 0.878 x lO^kg/m^, u = 3.94 x IQ-^rnVs. 




Figure 4: Critical current in dependence on 77 for Mg at 700°C, under the assumption of 
a current density of 10 kA/m^. 



Let us assume now a current density of 10 kA/m^ as it has been discussed 
as a typical achievable value for liquid metal batteries [1|. For Mg, the solid 
line in Fig. H] shows the increasing dependence of the critical current on the 
ratio Tj/ro. For ri/ro = the critical current is approximately 2kA (note 
that this value is independent of the radius). This means, in turn, that the 
maximum area vrr^ (for the assumed current density of lOkA/m^) is 0.2 m^. 
If we would like to increase the bearable current we have to increase rj/ro. In 
this case, however, the total current through the cylindrical gap also decreases 
as ~ (^o "''"?) 5 which is depicted by the various non-solid curves in Fig. |H At 
the crossing point between these curves with the solid curve one can identify 
the maximum current that would be possible for a given total area vrr^. For 
example, take the case of a rather large area vrr^ = 2m^. Evidently, for this 
case we cannot achieve an expected total current of 20 kA which would be a 
factor 10 higher than the critical current at ri/ro = 0. Yet we can get a total 
current of around 12.5 kA, if we choose r^/vo = 0.61. 

For sodium the situation is even more critical as shown in Fig. [51 Here, 
the current threshold at ri/ro = is approximately 0.86 kA, and even with 
an area vrr^ = 1 m^ we can only achieve a critical current of 6.1 kA if we go 
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Figure 5: Critical current in dependence on r\ for Na at 300°C, under the assumption of a 
current density of 10 kA/m^. 



to rj/ro = 0.62 

Note that the focus of this section was on identifying the threshold- 
increasing effect of hollow cylinders, without taking into account any (sig- 
nificant) axial current /a- It is clear, however, that the idea of the former 
section to avoid the TI by sending an opposite directed current through the 
middle remains valid also in the case of real-world liquid metals. 



5. Conclusion 

The main purpose of this paper was to point out that the Tl could repre- 
sent a serious obstacle for the upward scalability of liquid metal batteries. If 
it occurred in one of the liquid layers (most likely in the upper metal, having 
the highest conductivity and the lowest density) it would lead to a vigorous 
motion and very likely to a mixing of the layers with possibly "explosive" 
consequences. To give a first conservative estimate for the critical currents, 
we have used a simplified one-dimensional numerical model assuming an in- 
finite ratio of height to the radius of a battery with pre-supposed cylindrical 
geometry. For any realistic aspect ratio, this estimate has to be concretized 
by appropriate two-dimensional simulations [12], which would give larger val- 
ues of the critical currents. Another, although not a dramatic, modification 
of our results should be expected for the case that the shape of the battery 
is not a cylinder but, e.g., a cuboid. 
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We have also discussed two ways of how the TI can be avoided. One of 
them is the use of an inner non-conducting tube whose mere existence leads 
to a significant upward shift of the critical current. For a typical cell area 
in the order of 1 m^ this might indeed represent a viable technical means 
for excluding TI. A more radical way that can be expected to suppress the 
TI completely is just to return the current through the hollow cylinder in 
the middle. The same effect could be achieved by applying a parallel axial 
current which is at least a factor j^ as large as the current in the fluid. 

Many theoretical aspects of the TI in liquid metal batteries are yet to be 
discussed. First experiments are presently carried out to study the TI in a 
cylindrical liquid metal and to validate various possible ways how it can be 
avoided. 
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Appendix A. The mathematics of Tayler instabihty 

In this appendix we describe how to determine the critical current of 
the TI for the hypothetical case of a cylinder with infinite length. For this 
pur pose we use a simplified version of the procedure described in detail in 



10l |. We start with the Navier-Stokes equation for the velocity field V 



dV 11 

^— + (V-V)V = — VP + Z/AV + — (V xB) X B (A.l) 

Ot P jJ,Q 

and the induction equation for the magnetic field B 

rfVK 1 

^ = V X (V X B) + AB . (A.2) 

Here, p is the density of the fiuid, v its kinematic viscosity, a its electrical 
conductivity and pq the magnetic permeability of the free space. We also 
note that both fields have to be divergence-free (we assume incompressibility 
of the fluid): 

V-V = 0, V-B = 0. (A.3) 

Both V and B, as well as the pressure P, are now split into the basic state 
and a perturbation according to 

V = Vo + v P = Pq+P B = Bo + b. (A.4) 

In our particular case the basic state of the velocity is just Vq = 0, while the 
basic state of the magnetic field is given by 

zfl = (A.5) 

(A.6) 

(A.7) 

We employ now the usual normal mode analysis, by searching for solutions 
of the linearized equations for 

v(r, (f, z, t) = v(r) exp{i{kz + mip + ut)) (A. 8) 

p(r, ip, z, t) = p{r) exp(i(/c2; + rmp + ujt)) (A. 9) 

b(r, if, z, t) = b(r) exp{i{kz + rrup + ut)) (A. 10) 
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Br,o = 0, P^,o(r) = 


ar + /3/r, B 


with the constants a and /3 given by 


271 r'^ 1 — T]"^ 


^-t. 





This way, we end up with an eigenvalue system of ordinary differential equa- 
tions (in r): 

dvr Vr im , / . X 

-r^ + — + —V^ + lkv, = A.ll 

ar r r 



ar r V ^ 



+ i—X2 + ikXs + { k' + ^ ]u. 



-iB.ai^—B^flbr + 2B.ai^—^b^ = (A.12) 
~ {k^ + ^r\v^ + 2t—v. 



dr \ r"^ J ^ r"^ ^ 

+ —^ — {rB^^o)br + i}iai^—B^flb^-i—p = (A.13) 

-7^ + — -A;2 + — U, -iA;p + 2Hai2-5^,o6. = A.14 
dr r \ r"^ J r 

dbr br im 

— + - + - 
dr r r 

dbz i /, 9 m^\ , m ,^ m „ / . x 

-r-T{k^ + ^)br + ^X^-—B^,oUr = A.16 
dr k \ r'^ J kr kr 

dX^ / , 9 vp?\ ^ . Im , 



- + - + —b^ + ikb, = (A.15) 



dr \ r^ / r^ 



_r—{^^)vr + t-B^flV^ = (A.17) 
dr \ r J r 

where the abbreviations 

X, = ^ + ^, X3 = ^, X, = ^ + ^ (A.18) 

ar r dr dr r 

have been used. The governing parameter Hai here is the so-called Hartmann 
number, defined with the magnetic field at the inner radius according to 



Ha; = B^{r,)^r,{r,-r-;)^a/{pv) (A. 19) 

= ^Vri(ro-rOv/^7(H- (A.20) 

In the numerical implementation we will always use a (possibly very small) 

axial current. The boundary conditions, both at n and Tq, for the velocity 
are the usual no-slip conditions 

Vr = V^ = V, = Q (A.21) 
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As for the magnetic field, we assume insulating boundary conditions which 
can be written as 

K-^b, = (A.22) 

br + YJ]^J-^.^"^if'''0 - Irn+i{krO)] K = (A.23) 



at the inner boundary ri and 

h^ - 

kr, 



TTl 

b,-—h = (A.24) 



O 



br + -rnr^ -—K^ikro) - K^+iikr^)) ]h = (A.25) 

at the outer boundary r^. Note that /„ and K^ represent the modified Bessel 
functions. 

In order to determine the critical currents for a given ratio of rj and 
Ja//f we solve the above eigenvalue system by means of a shooting method 



adopted from [IJ]. Then we have to look for those wavenumbers k that lead 



to a minimum critical L 



References 



[1] D. Sadoway, G. Ceder, D. Bradwell, High-amperage energy storage de- 
vice and method, US Patent apphcation, US 2008/0044725 Al, 2008. 

[2] D. Bradwell, Technical and economic feasibility of a high-temperature 
self-assembling battery. Master Thesis, MIT, 2006. 

[3] B. Agruss, Regenerative Battery, US Patent 3,245,836, 1966. 

[4] T. Oshima, M. Kajita, A. Okuno, Development of Sodium-Sulfur Bat- 
teries, Int J. Appl. Ceram. Techn. 1 (2004) 269-276. 

[5] T. Sele, Instabilities of the metal surfaces in electrolytic aluminum re- 
duction cells. Met. Trans. B 8 (1977) 613-618. 

[6] P.A. Davidson, Overcoming instabilities in aluminium reduction cells: a 
route to cheaper aluminium. Mat. Sci. Technol. 16 (2000) 475-479. 



14 



[7] R.J. Tayler, Adiabatic stability of stars containing magnetic fields. 1 - 
Toroidal fields, Mon. Not. R. Astron. Soc. 161 (1973) 365-380. 

[8] Y.V. Vandakurov, Theory for the stability of a star with a toroidal 
magnetic field, Sov. Astron. 16 (1972) 265-272. 

[9] W.F. Bergerson, CB. Forest, G. Fiksel, D.A. Hannum, R. Kendrick, 
J.S. Sarff, S. Stambler, Onset and saturation of the kink instability in a 
current-carrying line-tied plasma, Phys. Rev. Lett. 96 (2006), 015004. 

[10] G. Riidiger, M. Schultz, D. Shalybkov, R. Hollerbach, Theory of 
current-driven instability experiments in magnetic Taylor-Couette 
flows, Phys. Rev. E 76 (2007) 056309. 

[11] G. Riidiger, M. Schultz, Tayler instability of toroidal magnetic fields in 
MHD Taylor-Couette flows, Astron. Nachr. 331 (2010), 121-129. 

[12] M. Gellert, G. Riidiger, Eddy diffusivity from hydromagnetic Taylor- 
Couette flow experiments, Phys. Rev. E 80 (2009), 046314. 

[13] B. Agruss, H.B. Karas, The thermally regenerative liquid metal concen- 
tration cell, in: C.E. Crouthamel, H.L. Recht (Eds.), Regenerative EMF 
Cells, ACS, Washington, 1967, 62-81. 

[14] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.F. Flannery, Numerical 
Recipes, Cambridge University Press, Cambridge, 1992. 



15 



